home *** CD-ROM | disk | FTP | other *** search
- /*
- ### product or two integer matrices ###
- */
-
- mprod(q,a,b,narow,nacol,nbcol)
- double *q, *a, *b;
- int narow,nacol,nbcol;
- {
- int i,j,m,kij,kim,kmj;
- for(i=0;i<narow;i++){
- for(j=0;j<nbcol;j++){
- kij=i*nbcol+j;
- *(q+kij) = 0;
- for(m=0;m<nacol;m++){
- kim=i*nacol+m;
- kmj=m*nbcol+j;
- *(q+kij) += *(a+kim) * *(b+kmj);
- }
- }
- }
- }
-